Response of Cloud Condensation Nuclei (> 50 nm) to changes in ion-nucleation 
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In experiments where ultraviolet light produces aerosols from trace amounts of ozone, sulphur 
dioxide, and water vapour, the number of additional small particles produced by ionization by 
gamma sources all grow up to diameters larger than 50 nm, appropriate for cloud condensation 
nuclei. This result contradicts both ion-free control experiments and also theoretical models that 
predict a decline in the response of larger particles due to an insufficiency of condensable gases 
(which leads to slower growth) and to larger losses by coagulation between the particles. This 
unpredicted experimental finding points to a process not included in current theoretical models, 
possibly an ion-induced formation of sulphuric acid in small clusters. 



The role of ionization in atmospheric processes has 
been a controversial matter since it was suggested fifty 
years ago[l[ [3], and found in the correlation between 
global cloud cover and the influx of galactic cosmic rays 
(GCR)t3j]- Subsequent studies have shown correlations 
between GCR variations and changes in aerosol counts 
and cloud properties in the atmosphere [JQ, but these 
are stiU disputed 

Fortunately, the issue can also be addressed in the 
laboratory. Experimental evidence for a microphysical 
mechanism was first reported in 2007^01] and further ex- 
periments have recently added to its credibility [ll|, 
These experiments initially showed that an increase in 
ionization leads to an increase in the formation of ultra- 
fine aerosols (w 3 nm), but in the real atmosphere such 
small particles have to grow by coagulation and intake of 
condensable gases to become cloud condensation nuclei 
(CCN) (> 50 nm) in order to have an effect on clouds 13, 
chapter 17]. 

Theoretical doubts about the likelihood of such par- 
ticle growth into CCN have arisen from consideration 
of (1) the competition between the additional ultra-fine 
aerosols for the limited supply of condensable gases lead- 
ing to a slower growth and (2) the larger losses of the ad- 
ditional particles during the longer growth-time to larger 
particles by coagulation and by other loss mechanisms. 
Indeed numerical studies using the current knowledge of 
aerosol dynamics predict that variations in the count of 
ultra-fine aerosols will lead only to an insignificant change 
in the count of CCN 14 , [l^ . It is even suggested that an 
increased production of ultra-fine particles as a result of 
GCR ionization leads to a reduction in the CCN count. 

In order to study the growth of aerosols to CCN sizes, 
measurements were performed in an 8 m^ reaction cham- 
ber (SKY2) made from electro-polished stainless steel 
shown schematically in FiglT] One side was fitted with a 
Teflon foil to allow ultraviolet light (253.7 nm) to illumi- 
nate the chamber, which was continuously flushed with 
dry purified air. Variable concentrations of water vapor 
(H2O), ozone (O3), and sulphur dioxide (SO2) could be 
added to the chamber, where the pressure was held a few 
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FIG. 1: Schematic diagram of the SKY2 experiment. 



Pa above atmospheric pressure, and the temperature at 
around 296 K. The UV- lamps initiated a photochemical 
reaction producing sulphuric acid (H2SO4). 

Ions were produced in the chamber by the naturally oc- 
curring GCR and by background radiation from radon, 
and the ionization could be enhanced with two Cs-137 
gamma sources (30 MBq), mounted on each side of the 
chamber. The total number of aerosols generated in the 
chamber were measured with a TSI Model 3775 Conden- 
sation Particle Counter (CPC) with a cutoff at 4 nm. 
A particle size measurement was done with an electro- 
static classifier (TSI model 3080) fitted with a nano-DMA 
(TSI model 3085) covering the range 3-65 nm and a CPC 
(TSI 3025A). Ultra-fine H2S04-water aerosols (ss 10 nm) 
could be generated with an Electrospray Aerosol Gener- 
ator (TSI model 3480). 

Concentrations of ozone were measured with Teledyne 
T400 analyzer and sulfur dioxide with a Thermo 43 CTL 
analyzer. The chamber was also equipped with instru- 
ments to measure temperature, differential and absolute 
pressure, humidity, and UV intensity. 

Estimates of the sulphuric acid concentrations were 
made by measuring the growth-rate of particle diameters 
above 3 nmil6!|. This method gave typical concentrations 
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FIG. 2: During a typical experimental run, number densities 
of aerosol particles of increasing diameter were observed as 
a function of time. Left panel:- 3-10 nm (black), 10-20 nm 
(purple), 20-30 nm (dark blue). Middle panel: 30-40 nm (light 
blue), 40-50 nm (green), 50-60 nm (yellow). Right panel: GO- 
GS nm (red). At ~ 19.2 days the gamma sources were opened 
to increase the ionization as described in the text, and an 
increase in aerosol density began immediately in the 3-10 nm 
curve (black). Subsequently the increase in number densities 
slowly propagated down to the larger aerosol sizes, as shown 
by the slanting arrow. Note that the number of particles in 
the first bin (3-10 nm) is relatively low because of a lower 
sensitivity of the instrument to the smallest aerosols. Black 
curves are an average over 67.5 minutes. 



in the range 1-10 ppt, i.e., ~ 10^ - 10^ molecules cm 

The experiments were run in a mode where steady 
state conditions of H2SO4 were achieved under contin- 
uous exposure to the UV- light. Typically the gas mix- 
ture in the reaction chamber consisted of 40-50 ppb O3, 
0.8-1.0 ppb SO2 and a relative humidity of 25%. 

Experiments were performed under various levels of 
ionization and UV intensity in the chamber. After chang- 
ing one of these parameters the aerosols were allowed to 
grow and conditions to settle to a new steady state for a 
period of 24-36 hours. For example the gamma sources 
were opened resulting in an increase in ionization from 
about 3 ions-pairs cm~'^s~"'^ to 60 ions-pairs cm^'^s"^. 
This resulted in an increase of about 20 % in the forma- 
tion of small aerosols. The parameters were then kept 
constant for a period of about 36 hours until the new 
steady state was achieved. 

Figure [2] shows such a run, where the number densi- 
ties of seven sizes of aerosols, from 6.5 nm to 64 nm, are 
plotted as a function of time. When the gamma sources 
were opened, at «19.2 days, an increase in aerosol density 
was seen to follow directly in the 3-10 nm curve (black) 
and then to propagate slowly through increasing aerosol 
diameters down to the larger aerosol sizes. After a new 
steady state was reached, the gamma sources were closed. 



and a small decrease in aerosol density could be moni- 
tored, again starting quickly with the small aerosols and 
propagating gradually to the larger aerosols (not shown 
in Fig. EJ. 

For each experimental run the density of particles be- 
fore (and after) an imposed ionization change was aver- 
aged over a period of 2.25 hours (prior to and after the 
change) and the mean and standard error of the mean 
was calculated. And finally the change in the response 
was averaged over five runs. The blue circles in Fig. 
show the relative response to changes in ion-nucleation 
as a function of particle size, averaged over the five runs. 
It is seen that the response is remarkably constant over 
the shown size range. 

It is of interest to contrast the above experiment to a 
situation without ionization and a constant H2SO4 pro- 
duction subject to an increase in ultra-fine aerosols ~ 6-8 
nm. The experimental procedure is first to reach steady 
state conditions using a constant UV intensity and trace 
gas concentrations as before, followed by a constant in- 
jection of Il2S04-water ultra-fine aerosols produced by 
the Electrospray Aerosol Generator. Figure |3Jd displays 
the response, averaged over 3.5 hours, as a function of 
aerosol sizes to the aerosol injection. It is seen that in 
this case the response is diminishing as a function of size 
in accordance with the theoretical expectations. 

The experimental results can be compared with nu- 
merical simulations of a general dynamics equation of 
aerosols. The evolution of the cluster distribution is given 
by [17 1 



dt 
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Nk + h-iNk-i - PkNk + SSk,ko (1) 



where Nk is the number density of clusters each contain- 
ing k sulphuric acid molecules, assuming that the equi- 
librium concentration of water molecules in each cluster 
is reached instantaneously 18j. K^j is the coagulation 
coefficient determined from Laakso et al. [l^, and can 
be used for all Knudsen numbers and hence from diame- 
ters of < 1 nm to > 1 microns. The radius of a cluster 
with i II2SO4 molecules and a number of water molecules 
depends on the humidity [isl. Chap. 10]. Particle loss to 
the chamber walls is approximated with the A/r^ term 
where 7 is determined experimentally to 7 =0.69 ± 0.05 
from the decay of particles in the chamber, and A — 6.2 
± 2.0 -10"^ nm'''~-'^s~^. S is the production of new criti- 
cal clusters with size given by fco II2SO4 molecules. The 
nucleation rate S is either a constant S ~ Sq or function 
of the II2SO4 concentration, e.g. S = a [112804]^, where 
a is a constant. The /3fe-term describes the condensation 
of II2SO4 molecules in the gas phase to the A;'th cluster 
and are found according to Laakso et al. 0|, with the 
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FIG. 3: Steady state response to a change in nucleation as a 
function of particle diameter, normalized to the particle num- 
ber before two types of perturbation, a) Ion-induced increase 
in nucleation. Blue circles are the experimental results av- 
eraged over five runs. The red curve is a typical result of 
a numerical simulation of the experimental situation using a 
standard numerical aerosol model. Notice that the expected 
response from the modeling decreases strongly with particle 
diameter in contrast with the experimental results. A much 
better agreement is seen with a numerical simulation in the 
black curve, where the concentration of sulphuric acid is held 
constant, b) Control experiment where the increase in par- 
ticle concentration is done by injection of Il2S04-water ultra 
fine aerosols (~ 6-8 nm) under constant UV intensity and 
trace gas concentrations, and no gamma source ionization. 
Notice that in this case the response (red circles) diminishes 
as aerosol size increases. Error bars are ±1-(t errors. 



value 1 of the mass accommodation coefficie nt 12Q] , and a 
mean free path from Lehtinen and Kulmala [21|. 

The equation governing the sulphuric acid concentra- 
tion is 



d[H2SOi] 
dt 



Ph.so, -{L + Xh,soJ[H2S04] (2) 



where PH2SO4 is the production of gaseous sulphuric acid. 
The second term L is the loss of H2SO4 molecules to the 
aerosols by condensation. The last term is the loss of 
H2SO4 molecules to the chamber walls, and is determined 
from extrapolating the size dependent aerosol losses to 
the size of a H2SO4 molecule to XH2SO4 = 7.2 ± 3.0 -lO""* 
s~^. The model is described in more detail in Bondo et al. 

The red curve in Fig. [3^ shows the result of a numeri- 
cal simulation of an increase in the nucleation similar to 
the experimental situation, where the production term of 
H2SO4 is kept constant, but the concentration of H2SO4 



can vary. The response goes slowly to zero with increas- 
ing size of the clusters, due to a smaller concentration of 
H2SO4. This simulation of the aerosol dynamics is con- 
sistent with a number of recent simulations which show 
very small responses at CCN sizes to a change in the nu- 
cleation rate. In their MODGIL simulation Pierce and 
Adams found responses of 0.004% in CCN (at 0.2% su- 
persaturation) to a 4-fold increase in new particle forma- 
tion, and in another simulation (lONLIMIT) they found 
a 0.08% change in CCN (again at 0.2% supersaturation) 
to a 24% increase in nucleation 

As the expectation shown in the red curve in Fig. [3^ is 
contradicted by the experimental results (blue circles) an 
obvious question is whether the ionization by gamma rays 
may produce sufficient H2SO4 in the gas phase to replace 
the expected loss of « 3-7% from the additional parti- 
cles (estimated from numerical simulations). Each ion- 
pair will on average 23. 2^ produce two OH molecules. 



Therefore with an ionization of 60 ion-pair cm~'^s~^ the 
production will be 120 molecules cm~^s~^. From the ex- 
perimentally estimated losses and growth rates the pro- 
duction of H2SO4 from the photolysis is 3.5T0'* molecules 
cm~'^s^^. If every OH molecule becomes an H2SO4 
molecule its production will be 0.3% of the photolysis, i.e 
10 times lower than estimated loss of H2SO4. From the 
experiment it is known that only a minor fraction of the 
OH are consumed in the pathway that leads to H2SO4. 
It is therefore safe to conclude that the production of 
H2SO4 by this path is at least an order of magnitude too 
small to explain the observed. Another main chemical 
species produced in moist air (N2-O2-H2O) is nitric acid 
(HNO3) that potentially could help the condensational 
growth of the aerosols. Experimentally it is found that 
each ion-pair produces 0.4 HNO3 molecules [25|. which in 
the chamber leads to a production of 24 cm~'^s~^ HNO3 
molecules. This production is only 0.07% of the photol- 
ysis production of H2SO4, and therefore nearly 40 times 
to small too explain the experimental results. 

In a second simulation, illustrated in the black curve 
in Fig. [3^1, the concentration of H2SO4 is artificially held 
constant. In this case the response of larger particles 
to the additional nucleated particles is not going to zero 
and the match to the experimental results is much bet- 
ter. How, then, is the growth of the particles sustained? 
The indication from the second numerical simulation is 
that effectively there is no decrease in the concentration 
of condensable gases, even though the UV photolysis of 
H2SO4 is held constant throughout the duration of the 
experiment. But the additional ion-nucleated particles 
should effectively decrease the H2SO4 concentration with 
w 3-7 %. 

A possible explanation could be that the charged clus- 
ters are producing additional H2SO4 molecules from re- 
actions involving negative ion chemistry of O3, SO2 and 
H2O, where a negative ion can be reused in a catalytic 
production of several H2SO4. Such reactions were first 



4 



suggested in Svensmark et al. [lO|, and also in a recent 
experiment [2^ looking at isotope fractionation of sulphur 
from either UV or from ion- induced generation of H2SO4, 
where the sulphur isotope fractionation was used to dis- 
tinguish the different pathways leading to H2SO4. It was 
found in the presence of ionization alone, that for each 
ion-pair 27.8-10^ H2SO4 molecules were produced, using 
extreme gas mixing ratios, i.e. 0.01 % SO2, 400 pbb 
O3, 40 % RH H2O and 1000 ions-pairs cm-^g-i. Scal- 
ing to the present experiment gives a production of 38.9 
H2SO4 molecules pr. ion-pair, and so a total production 
of w 2.3-10^ H2SO4 molecules cm'^^s"^. This amounts 
to an increase of « 7 % in the production of H2SO4 
molecules which is sufRcient to compensate for the nu- 
merically determined w 3-7 % decrease in H2SO4 con- 
centration caused by the increase in number of aerosols. 
Recent ab inito calculations have confirmed the first steps 
of such reactions in small clusters pTj - fiof . 

So in conclusion it has been shown that an increase in 
ion-induced nucleation survives as the clusters grow into 
CCN sizes in direct contrast to the present neutral ex- 
periment and current theoretical expectations. It is pro- 
posed that an ion-mechanism exists which provides a sec- 
ond significant pathway for making additional H2SO4, as 
a possible explanation of the present experimental find- 
ings. Irrespective of the detailed mechanism leading to 
the results presented here they provide a possibly impor- 
tant missing piece of the puzzle as to why responses in 
aerosol to variations in ionization have been seen in cloud 
properties. 
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